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Abstract 


In the last years many countries in the world have opened their power 
markets to allow competition between power producers At the first 
look it seems that the renewable energies do not have a chance in 
those deregulated markets because they are more expensive than the 
conventional one and some of them are having problems because of 
fluctuations in the production In this thesis the fluctuation problem is 
addressed and forecast models for the production of renewable energies 
mainly wind and photovoltaic are developed This will help the trading 
of these lenewable energies m deregulated powei market For 
forecasting the production of renewable energies fuzzy neural network 
and regression models are utilized and tested with the actual measured 
data 
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1. Introduction 


In the last years many countries in the world have opened their power markets, to 
allow competition between power producers At the first look, it seems that the 
renewable energies do not have a chance in those deregulated markets, because 
they are more expensive than the conventional one and some of them are having 
problems because of fluctuations in the production 

The mam advantage of the renewable energies compared with the conventional 
energies is that they are not causing any appreciable environmental disturbances, 
like C02 increase and greenhouse effect Because of these advantages the pro 
duction of renewable energies is subsidized in some of the industrial countries, 
which already are deregulated The purpose of these subventions, is to make the 
renewable energies competitive at the power market 

Although the renewable energies were made competitive, then position in the 
deregulated market would still not be simple, because of fluctuations in the pro 
duction This applies specially to wind and photovoltaic energy To make the 
trading of those ’’fluctuating renewable energies” possible, forecasts for their pro- 
duction is obviously of great importance With such a forecast, the buyers of 
renewable energies, have the chance to compensate the coming fluctuations m 
the renewable production through trading at some conventional power markets, 
for example the spot market 

In this work some forecast models for renewable energies will be designed, im 
plemented and compered The models will be limited to wind and photovoltaic 
energy that is the most fluctuating ones For the wind power forecast a model 
based on fuzzy logic will be compered with a regression analysis model For the 
photovoltaic power forecast a neural network model will be established and also 
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compered with regression analysis model 

In next chapter the trading of renewable energy in deregulated market and the 
need for power forecasts will be discussed further Chapter 3 and 4 handle the 
mam theme of this work, the wind and the photovoltaic power forecasting mod 
els Chapter 5 and 6 deal with the basics for multiple regression and neuronal 
networks, that is theory that is used for the wind and the photovoltaic models 
The conclusions of the work are finally in chapter 7 


2 The Trading of Renewable 
Energy in a Deregulated Market 


In this chapter the problematic concerning the trading of renewable energy m a 
deregulated market will be discussed The two mam problems that a renewable 
energy maiketei has to deal with aie 

• The expensiveness of the renewable energy 

• The fluctuations of the renewable pioduction 

In practice there already exist some techniques to deal with the first problem, 
for example the renewable order in England and Wales [17] and the green label 
approach in Holland [22] Those techniques have in common that their aim is to 
make the renewable energy competitive at the power market through some kind of 
furtherance The second problem, the fluctuations of the renewable production, 
is harder to handle A power maiketer, who is purchasing or selling renewable 
energy, has to take those fluctuations into account 

Heie the focus is on the second pioblem Foi that purpose a simple model of a 
deregulated powei maiket will be set up and the situation of a power marketer 
who IS trading with renewable energy, will be analysed It is admitted that the 
renewable production is promoted in some way, so it is competitive It will be 
shown that the marketer will definitely need a forecast for the production of the 
renewable energy 
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Figure 2 1 A model of a deregulated power market 


2 1 A Model for a Deregulated Power Market 

Figuie 2 1 shows a simple model of a deregulated power market The market is 
strictly divided into 4 sectors Production, transport, marketing and distribution 
The power production is in the hands of independent power producers [IPP’s) 
They are classified as renewable and conventional A grid operator is responsible 
for the transport of power from the IPP's to the customers The grid operator 
does not produce any power at his own He will buy regulation power directly 
from some conventional power producers The spot market is an organised day 
ahead market for electric energy [20] The participants at the spot market give 
their offer how much they want to buy/sell and for what price the day before 
the activation of trade A power marketer buys power from the IPP's and sell 
it to his customers The customer do not handle directly with the IPP' s It 
IS assumed that the analytical method [5] is used, that means that the power 
maiketci has to pay the difference between the consumption of his customer and 
the estimated demand A power marketer has 4 ways to buy energy for his 
customers, as indicated in figure 2 2 Those 4 ways are 

1 Buying energy through bilateral trading from conventional IPP 

2 Buying energy through bilateral trading from renewable IPP 

3 Buying energy through the spot market 

4 Buying regulation power from the grid operator 
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Figure 2 2 The possible ways foi a power marketer to buy eneigy for his cus 
tomeis 

If the power marketers haven’t order enough power through bilateral trading or 
through the spot market to fulfil the consumption of their customers , the grid 
operatoi has to buy regulation power at the regulation market to balance out the 
difference betweeir consumption and production The power marketers have to 
pay for the regulation power 

2 2 Trading of Renewable Energy 

Now the situation of a power marketer, who is trading with renewable energy, 
will be considered It is assumed, that the marketer is having contract with 
one or more IPP' sren and that he is buying all the power that they can 
possibly produce Moreover it is admitted that the power marketer is buying 
some constant amount of power from one or more IPP' scoN He wants to buy 
the rest, that he needs to fulfil his customers consumption, from the spot market 
Figure 2 3 shows this for a one particular day The power marketer has to decide 
the day before this particular day, how much he wants to order from the spot 
market for each period of the following day Here each day is divided into 24 
periods, that is, each period is 1 hour What the power marketer needs from the 
spot market in some particular period is, as indicated in figure 2 3, the estimated 
consumption minus what he will get from the IP P's It’s clearly important that 
the marketer has some good forecasts for the consumption of his customers and 
the production of the IPP^^^s for the following day Here (in Figure 2 3) is 
admitted that he has such forecasts, with constant values during each period 
Figure 2 4 shows the difference between the estimated and real consumption for 
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Figure 2.4.: Estimated and real consumption 
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this one particular day The grid opeiator balance out the differences (positive 
or negative) and the power marketers have to pay for it afterwards So good 
estimation for the consumption is not just helpful, when ordering from the spot 
market, but also extremely important because of the regulation costs 

2 3 Needed Forecasts for a Power Marketer 

As a result of this chapter the following can be stated A power marketer who 
IS responsible for selling power from some IPP'sren, and who wants to com 
pensate the fluctuations of the renewable production with purchase at the spot 
market, needs the following forecasts for the coming day 

• The estimated consumption for each period 

• The estimated production of the IPP' Sren for each day 

In this report the second problem, that is the estimated production of the IPP'sren 
for each day, will be considered This work is concentrated on wind and photo 
valtiac power forecasts, because those are the renewable energies were one can 
expect the most fluctuations 



3 Wind Power Forecasting 


3 1 introduction 

Wind energy converters (WEC) integrated in electrical power systems may cause 
several pioblems, including power quality, stability and especially power dispatch 
ing These pioblems increase as the penetration of wind energy increases The 
power generated by wind turbines changes rapidly because of the continuous 
fluctuation of wind speed and direction As explained m previous chapter it is 
extremely important, for the trading of wind energy m deregulated power market, 
to have good forecast for the wind power production at any specified time for the 
following day 

In the literature, various models, have been developed for wind speed or power 
prediction The time scale of these models ranges from some seconds, when the 
objective is wind turbines control, to some minutes, or even hours, when the ob 
jective IS economic dispatch and powersystem planning The time series analysis 
is common model for small time scales In higher time scales, of some hours or 
more, models incorporating meteorological information are usually used The 
time series models found in the liteiatuie do not provide a significant improve- 
ment over the method of persistence As a consequence, the latest one, being 
extremely simple (the wind in all future time steps will be equal to the wind 
now), IS almost exclusively used m the applications [15] In [6] ARMA (auto 
regressive moving average) models for wind speed prediction achieve an improve- 
ment with respect to the peisistance model In [4] multilayer feed-forward neural 
networks and radial basis functions are applied for wind power prediction In [9] 
various models like ARMA and bilinear ones are developed for wind power In 
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[1] an artificial neural netwoik based model which uses local and spatial relations 
of the wind speed is proposed foi estimating wind power 

In this chapter piobability and fuzzy ba^ed model is developed for short term 
piediction of wind power This model uses wind velocity and its regularity to 
forecast the wind power The basic idea is developing a possibility cuive for wind 
powei generation at a forecasted wind speed and apply some fuzzy operations, 
like fuzzification, composition and defuzzification, to get an approximate wind 
power Wind power possibility curve is obtained by manipulating the historical 
data At last regression analysis is also used to make comparison for the proposed 
model 

This chapter is divided into 5 sections In section 3 2 the basics of the forecasting 
methodologies aie explained Section 3 3 explains the detailed description of the 
proposed model where as section 3 4 discusses about the regression model In 
section 3 5 results are presented for the proposed fuzzy model as well as regres 
Sion model and compared both In section 3 6 conclusions regarding the proposed 
model with lespect to the regression model are discussed 

3 2 The Basics of Wind Power Forecasting 

The power generated by a wind turbine is chiefly influenced by wind velocity 
and direction [26] However, many other factors such as the turbines immediate 
histoiy of power generation, the blade position at the time of measurements [21] 
dynamic load distiibution among paiallel tui bines, alternation of system load 
flow, dynamic perfoimance of a generator, and wind dynamics also influence gen 
erated powei, so that even when the average wind speed within 10 minutes is the 
same, the power generated can be very different [24] 

The diiection of wind also influence power generation However, compared with 
wind velocity, wind direction has less influence on power output [24] because most 
of the turbines are built to turn into the wind when operating Generally, at the 
same level of wind speed, there is no great difference of the power generation for 
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different wind direction 

Theie aie two main approaches that may be followed in order to generate wind 
power forecast 

1 to develop an explicit prediction model for wind power, in which it will 
be possible to consider wind speed, wind direction etc as explanatory 
variables 

2 to develop a prediction model for wind speed and a second model for the 
transformation of wind speed to power 

In this project the first methodology is used for forecasting wind power in which 
wind speed and wind speed regularity are considered as input parameters The 
detailed model is discussed in next section Here are presented Naive predictors 
foi the wind power forecasting in which only data manipulation is involved 

3 2 1 The Naive Predictors 

When an advanced forecasting model is not available, forecast may be obtained 
with a minimal effort and data manipulation and can be based solely on the most 
recent information available Such forecasts are referred to as Naive forecasts 
One such method (Persistent or Naive 1 method), is to use the most recent datum 
available P(t), as forecast P{t k/t) for each one of the future time steps that 
IS P(t4-k/t) = P(t) for fe = 1,2, n A slightly more sophisticated method 
would be to use the average of m past values (Naive m) as forecast 
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Figure 3 1 Fuzzy based model 

3 3 The Fuzzy Based Model 

This model can be divided into two parts One is the probability part and 
the other is the fuzzy part as shown in figure 3 1 In the probability part a 
probability curve of the wind power generation at a forecasted wind velocity will 
be developed by using historical data In the fuzzy pait the probability curve will 
be manipulated through some fuzzy opeiations according to information about 
the regularity 

3 3 1 The Probability Part 

In this section developing probability curve of the wind power generation at a 
forecasted wind speed is processed This is done by using the historical measure- 
ments that are taken for wind speed and power generation As shown in figure 
3 1 the probability part is having two inputs One is historical data and the 
other one is forecasted wind speed Output of this block is the probability curve 
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of the wind powei In the data manipulation block, some statistical analysis of 
histoiical data with lespect to the forecasted wind speed will be carried Now 
the data manipulation process will be explained It is admitted that x is wind 
speed and y the corresponding wind power 

Consider theie are N pairs of measurements (®i, j/i)? (312? 2/2)5 li^ruVn) 
between two variables x and y are in historical data The first step in data ma 
mpulation block is analyzation of the relation ship between x and y As indicated 
m figure 3 2 , there may be several values of wind power for one value of wind 
speed As mentioned before that the output of the data manipulation process 
is a probability function for the wind power This probability function indicates 
the different values of wind power that one could expect for the forecasted wind 
speed according to the historical data Now the composition of this probability 
function will be be discussed 

Foi the composition of the probability curve only the points in the histoiical 
data, which lay between Xq —■ A® and Xq + Aa;, will be used (see figure 3 2 ) 
Those points are (aji, j/i), (01252/2)5 5 (®M5 2 /m)) where M is the number of 

points laying m this range Here Xo is the forecasted wind speed and Ax is the 
clustering range around the Xo 

Selection of the size of Ax depends on the density of the data points around Xp 
If the density is high. Ax should be small and vice versa This is reflected m the 
below two rules for Ax 


1 Ax should not exceed certain range that is Axmax-, and 

2 The number of points (Af) should not exceed certain value that is J\d.pnax 


Axmax and Mmax are constants, which the user can vary depending on his 
experience with the model 

After getting the points, it is better to do normalization to the above points, by 
using below equation 


^t{norm) 


Xi^origtnal) 

max(xi,X2, 


X 100 


(3 2) 
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Figure 3 2 Wind speed versus wind power graph 


Here a3t(norm) is the noimalized value and Xt(^ortgtnai) is the original value This 
normalization is for making the computation easy 

The probability function is nothing but the comparison among the points, so 
here comparison between the above points is done with it’s mean, which can be 
calculated by using below equation 


«Cc — 
Vc - 


Xi(narm) "F X2{norm) ”1” 


M 

2/l(uorTTi) "F 2/2(Tiorm) “F 


M 


“F XM(norTn) 


”F Dm (norm) 


(3 3) 
(3 4) 


Here Xc and yc represents the mean point among the M points 
The distances of these points from its mean is the one way to compare them 
Here the distances of each point from its mean is calculated and is denoted with 
di in below equation 


d^ 



(3 5) 


X%(norm))^ ”F (?/c Viinorm))^ 
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Figure 3 3 Probability of wind power at forecasted wind speed 


One can observe from the figure 3 2 that for the different wind speeds, there may 
be the same wind power It means that, in the above points there may be same 
value of y for different x values Here the process is to find the probability of 
wind power that means probability of getting each different wind power for a 
forecasted wind speed If there is some points which having same wind power, it 
could be better to find the average distance of those points, which in turn leflects 
the average vaiiation of this wind power to the mean wind power Consider there 
aie L diffeient wind power values, m those Af points The average distance of 
each different wind power can be obtained by below formula 

cii = ^ dr, (3 6) 

^ 1=1 


where di is the average distance of 1 th wind power and K is the same value of 
wind power points Now the probability of each wind power (i/z), that is p{yi) 
can find by using below equation 

100 - di 


PiVi) 


(3 7) 


100 
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After finding all the probability values for different wind power values, one can 
get the piobabihty curve of the wind power generation at a forecasted wind speed 
as shown in figure 3 3 

3 3 2 The Fuzzy Part 

This fuzzy part is having two inputs, one is the above probability function of 
the wind power at a forecasted wind speed and wind speed regularity informa 
tion The output of this fuzzy part is the estimated wind power for a forecasted 
wind speed Wind speed regularity information tells about the accuracy of the 
forecasted wind speed In this fuzzy part fuzzy operations like fuzzification of 
input variables, fuzzy composition and defuzzification is carried, which will be 
explained now 

3 3 2 1 The Fuzzification Process 

Fuzzification is carried for the input variables, that is "probability function of 
wind power generation” and "wind speed regularity information” Fuzzifica 
tion means converting the input variables into fuzzy representative membership 
function format Heie different fuzzification processes are applied for the input 
variables Next in this section fuzzification process for each input variable is ex- 
plained 

Fuzzification of wind speed regularity 

Regulaiity itself is a fuzzy term, so one can take wind speed regularity infoi 
mation as a fuzzy representation of natural language qualifiers "Very Regular”, 
"Regular”, and "Little Regular” displayed in figure 3 4 It can be represented in 
a functional form as shown below [18] 

fix) = 1 - (1 - (3 8) 

The values of </> are presented in Table 3 1 

Fuzzification of wind power probability curve 

Fuzzification of wind power probability curve is done by using normalization tech 
niques Here normalization technique that is applied for getting fuzzy wind power 
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Figure 3 4 Qualifiers for wind regularity 



Statements 

1 5 

Very Regular 

1 

Regular 

0 696 

Very Little Regular 


Table 3 1 Values of (j) foi wind speed regularity 

generation is, dividing all the probability values with its maximum probability 
value It makes the probability curve of the wind power into possibility curve of 
wind power generation as shown in figure 3 5 

3 3 2 2 The Fuzzy Composition 

Now there are two membership functions, one is fuzzy wind power generation 
and other one is wind speed regularity membership function So one can obtain 
modified fuzzy power generation by applying fuzzy composition rule of inference 
[19] In [18] traditional modus ponen rules are used to modify the fuzzy power 
generation 

Modus ponens deduction is a very common inference scheme used in forward 
chaining rule-based expert systems [19] It is an operation whose task is to find 
the truth value of a consequent m a production rule, given the truth value of 
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Figure 3 5 Fuzzy membership function of the wind power generation 

the antecedent in the rule Modus ponens deduction concludes that given two 
propositions P and P Q both of which are true, then the truth of the simple 
proposition Q is automatically inferred by below equation 

{PA{P-^Q))-^Q (3 9) 

In [19] different types of fuzzy composition rules are discussed in which ’’max 
min” composition is mostly used 

The final fuzzy power generation can be get by using ’’max min” composition 
which IS presented below 

tiPoM{x,y) = s\ipmin(np{x,y),fXM{y z)) (3 10) 

ver 

where x is fuzzy value for power 
y is membership value for the basic models, and 
z IS modified possibility values 

When a qualifier such as ’’little regular” is applied (curve LR of figure 3 4) , low 
values of possibility rise, and the resulting distribution is more ’’wide” , meaning 
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Wind PowerCp u) 

Figure 3 6 Modified wind power membership function 


that we are less certain about the true value than in the base case This can be 
shown graphically as shown in figure 3 6 If the qualification is regular then the 
original function remains unchanged 

3 3 2 3 The Defuzzification Process 

After getting the modified fuzzy power generation, one can get the crisp value 
of powei by using the fuzzy defuzzification techniques In liteiature [16] many 
defuzzification techniques aie discussed, but the ’’centroid method” is the most 
popular one 

The ’’centroid method” will now be explained, with the help of figure 3 7 In 
discrete case where x = [aii, 0:25 j ®i] and corresponding membership values 
fjL{x) = [fjLi{x)^rn '^ 2 {x), , /xi(aj)], then by using centroid method the crisp 

value of X is 

(3 11) 


X* — ^ 1=1 

eLi 



Chapter 3 Wind Power Forecasting 


Page 19 



Figure 3 7 Centroid defuzzification method 


By employing ’’Centroid” defuzzification method, for the modified wind power 
membership function one can get the estimated wind power for the forecasted 
wind speed 

3 4 The Regression Model 

Regression based prediction model also can be used for wind power forecasting 
Regression analysis involves basically statistical analyzation of the historical data 
One can write the regression equation for wind power forecasting as shown below 

Wzndpower = a + b{\Vindspeed) (3 12) 

where a and b are the regression equation coefficients, which are to be determined 
by using historical data The detailed description of the regression technique is 
explained in chapter 5 

3 5 Results and Comparisons 

The data that is taken to test the proposed model is from the ’’Norden Ostermarch” 
wind power plant in Germany They measures the data with a 5 minutes inter- 
val in each day and they measures mainly wind velocities and wind power The 
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measurement of wind speeds are taken at a hight of 10m from the sea level The 
rated powei of this plant is SOOkW, and rotor diameter is 37m which is located 
40m high fiom the sea level 

The data taken from the Norden Ostermarsch, for testing the above both models 
IS 5 minutes interval data of wind velocities and wind power generation from 11th 
Octobei to 14th November 1999 In this data the first one week data that is 11th 
October to 17th October, is used as a historical data and two other days (8th 
Octobei and 14th November) are selected for testing the model The analysis 
of the historical data is shown in figure 3 8 In this figure the first one is wind 
velocity versus wind power graph and other two are time versus wind speed and 
time veisus wind power graphs The same analysis is also carried for the testing 
data files and are shown m Appendix, figure A 1 and figure A 2 

3 5 1 Forecasting Results by Using the Fuzzy Based Model 

The above data is used to test the proposed fuzzy model and the results are 
shown in figure 3 9 and figure 3 10 corresponding to the 18th October and 14th 
Novembei data Although there is no information about regularity, heie the 
testing process is carried for all the three cases of the regularity information that 
IS very regular, regular and little regular The results with regularity information 
such as veiy regular and little regular are shown in Appendix, figure A 3 and 
figure A 4, foi 14th November data The figures figure 3 9 and figure 3 10, which 
aie aheady mentioned are corresponding to the ’’legulai” regularity information 
In Table 3 2 estimation error of the proposed model is presented with the actual 
and estimated value of wind power Error can be calculated by using below 
equation 


Wind Speed(m/s) 








PoweitkW) Power{kW) 
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Figure 3 9 18th October Fuzzy based model results 



Figure 3 10 14th November Fuzzy based model results 
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S No 

Wind Speed 
m/s 

Wind Speed 
Regularity 

Actual Power 
kW 

Estimated Power 
kW 

Error 

1 

32 

VR 

48 6 

45 75000 

0 0586419 

2 

47 

LR 

63 5 

111 417988 

0 754613 

3 

2 7 

R 

34 4 

36 226521 

0 0530965 

4 

5 5 

LR 

69 3 

150 363337 

1 169745 

5 

6 1 

R 

104 9 

175 042783 

0 668663 

6 

46 

VR 

90 2 

122 5 

-0 3580931 

7 

6 3 

VR 

58 5 

186 6 

2 1897435 

8 

38 

LR 

52 7 

66 506545 

0 261983 

9 

24 

VR 

20 9 

36 1 

0 727272 

10 

49 

LR 

124 2 

134 715801 

0 084668 


Table 3 2 Wind power forecasting by using fuzzy model 

3 5 2 Forecasting Results by Using the Regression Model 

Linear regression analysis is also used in the similar manner as the fuzzy model 
for forecasting the wind power The data from the 11th October to 17th October 
is used for the regression analysis and the other two data files are used for testing 
the regression technique Here wind speed is the only input variable taken into 
account for estimating the wind power 

Testing results of 18th October and 14th November data files are shown in figure 
3 11 and figure 3 12 Table 3 3 shows the results of estimated power for specific 
wind speeds with the Error The complete statistical analysis of the data bank 
with regression analysis is shown in Table 3 4 In that table the regression coef 
ficients a and b are presented with the multiple regression correlation coefficient 
R The standard error of b and the hypothesis test results for b and R that is 
Fb and Fr are also presented m that table Fb and Fr are the F test results, 
which IS a common test in statistics 


3 5 3 Comparisons 

By observing the fuzzy based model results with the regression based model, it is 
found that the regression technique is better than the fuzzy model Better results 










Power(kVV) Pawer(kW) 
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Figure 3 11 18th October multiple regression results 



Figure 3 12 14th November multiple regression results 
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S No 

Wind Speed 
m/s 

Actual Power 
kW 

Estimated Power 
kW 

Erroi 


3 2 

48 6 

22 508814 

0 53685564 


4 7 

63 5 

63 041298 

0 007223655 


2 7 

34 4 

08 997984 

0 73843074 

4 

55 

69 3 

84 658640 

0 22162531 

5 

6 1 

104 9 

100 87162 

0 03840211 

6 

46 

90 2 

60 339140 

0 33105162 

7 

63 

58 5 

106 27596 

0 8166831 

8 

38 

52 7 

38 721806 

0 26524088 

9 

24 

20 9 

0 8914852 

0 9573452 

10 

49 

124 2 

68 445640 

0 44890788 


Table 3 3 Wind power forecasting using multiple regression 


Parameter 

Regression 
Coefficient (6) 

Standard 
Error of b 

Hypothesis 
test Ft, 

Wind Speed 

] 0202 

0 0134 

5827 2 

intercept a 

Correlation Coefficient R 

-0 1569 

0 86206 


Hypothesis test Fr 

5827 2 



Table 3 4 Statistical analysis of wind data 


obtained from the regression model is can be explained by the statistical Table 
3 4 In that the correlation coefficient value R is very close to 1, signifies that 
the regression equation is best fit for the historical data The standard error of b 
IS also very less The F test results for both regression coefficient and coirelation 
coefficient shows that the both aie passed the null hypothesis test The detailed 
explanation regarding this tests are carried in chapter 5 

One mam advantage from the fuzzy model is, it can provide the possibility curve 
of the wind powei foi the foiecasted wind speed, which will provide the better 
sight about the wind power at that forecasted wind speed The results from 
the regression model are promissing, but the error is sometimes higher than 
accesptable value For more better results, one can modify the regression model 
by changing the regression equation from linear to the nonlinear as given below 

Wtndpower = a + 6 X Wtndspeed^ (3 14) 




















Power(kW) 
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This may give improved results, because it is found in literature that the lela 
tionship between wind powei and wind speed is non linear nature as mentioned 
above For the proposed fuzzy model one have to go further study on this model 
to make modifications for getting better results 

3 6 Conclusions 

In this chapter wind power forecasting models were discussed Two models were 
constructed for wind power forecasting, a fuzzy based model and a regression 
based model In fuzzy model, historical data manipulation is carried in the prob 
ability part, and in fuzzy part, fuzzy operations like fuzzification, composition 
and defuzzification were used Both the models were tested with the actual wind 
power plant data From the results it is concluded that the regression based 
model gave better results than the fuzzy model 


4 Photovoltaic Power Forecasting 


4 1 Introduction 

’ Clean electricity with no moving parts” is the powerful concept of photovoltaic 
powei systems Now a days photovoltaic systems are rapidly expanding and in 
creasing roles in electric power technologies, providing more secure power sources 
[25] The photovoltaic systems are, by nature, non linear power sources, that 
need accurate estimation of the maximum power generation for the operation as 
well as trading in deregulated powei market 

In [12] estimation of photovoltaic power by using neural networks is developed, 
which uses environmental factors such as irradiation, temperature and wind ve 
locity for the prediction of photovoltaic power generation In [13] the effect of 
clouds on photovoltaic power generation is studied 

The accurate prediction of the maximum power from the Photovoltaic power sys 
terns IS inevitable, because of high uncertainty of solar radiation The power from 
the Photovoltaic systems depends on environmental factors, mainly the iiradia 
tion and the cell temperature It is very hard to measure the cell temperature, 
but related factors, like the temperature at the solar module and the wind ve 
locity can be used instead In this study these factors along with the irradiation 
will be used foi the forecast 

For the prediction of maximum power generation, a three layer artificial neuial 
network is proposed The input signals are the irradiation, the temperature, and 
the wind velocity The output signal is the predicted maximum power from the 
photovoltaic system For the training data sets of the proposed neural network, 
two data sets, one for summer and other for winter are selected The neural 
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network has been trained by using one data set corresponding to one season, 
therefore, two types of the neural networks has been obtained from those data 
sets The accuracy of the prediction has been evaluated by using the actual data 
measured on the photovoltaic system The proposed model is also compared with 
conventional multiple regression analysis 

This chapter is divided into five sections in which section 4 2 discusses the basics 
of the photovoltaic power generation and solar cell characteristics Section 4 3 
deals with the forecasting methodologies mainly basics of the forecasting tech 
niques, proposed neural network model and regression model In section 4 5 
results and comparisons are discussed and in section 4 6 conclusion about the 
proposed model with its accuracy is mentioned 

4 2 Fundamental Concepts of the Photovoltaic 
Modules 

Before the description of the forecasting approach the fundamental concepts of 
photovoltaic power generation and solar cell operation will be discussed 

4 2 1 Photovoltaic Power Generation 

The generation of electricity from solar energy can be achieved through two major 
technology alternatives One uses the light from the sun to generate electricity 
directly (photovoltaic technologies), and the other uses the heat from the sun to 
mciease the tempeiatuie of a working fluid, which-in turn can be used to gener 
ate electricity (solai theimal technologies) Each of these major alternatives can, 
in turn, be subdivided into variants of the majoi technology Photovoltaic tech 
nologies fall into crystalline, multi crystalline, thin film or concentrator variants 
while the solar thermal technologies fall into trough, power tower, dish engine 
and thermal electric variants 

Photovoltaic power generation process is shown in figure 4 1 Generally speaking, 
photovoltaic solai cells use a semiconductor material that is exposed to sunlight 
The energy of the incident light displaces electrons from their normal atomic or 
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Figure 4 1 Photovoltaic power generation process 

bits and an electrode grid structure on the surface of the semiconductor, collects 
these electrons and makes them available for use in an external circuit This is 
very similar to the way that the chemical reaction and the electrodes in a dry 
battery cell make electrons available for external use The terms crystalline, thin 
film and concentrator desciibe the manner in which the semi conducting material 
IS processed and optimized as a photovoltaic cell Crystalline cells are fabricated 
from ingot of the semiconductor material, usually silicon, that are cut into rel 
atively thin slices, processed to optimize the electron collection efficiency and 
laminated into a protective enclosure Thin film cells are extremely thin layeis of 
semi-conducting material that are evaporated onto a substrate, and concentrat 
ing cells use a plastic lens to concentrate sunlight from a large area onto a much 
smaller area of ciystalline semi conducting material All types have their merits 
and problems and are described m detail in [27] Now after getting idea about 
the photovoltaic power generation, one have to understand about the solar cell 
characteristics which is described in next section 
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Rsh 


Figure 4 2 Equivalent circuit for a solar cell 


4 2 2 Solar Cell Characteristics 

A solar cell is simply a diode of large area forward bias with a photo voltage The 
photo voltage is created from the dissociation of electron hole pans created by 
incident photons within the built in field of the junction or diode The electrical 
behavior of the cell is describable by the simple equivalent circuit shown in figure 
4 2 Here the current- volt age characteristics for a solar cell are give by [14], 


I 

where Id 


J-L — J-D — 

Rsh 

^ ( fq(V + R.J)\ __ 

C’’*’ [- nkT ) - " 

^ , f fq{V + R,iy 

C"*' ( nkf • 


- 1 


V + R,I 

Rah 


(4 1 ) 
(4 2 ) 
(4 3 ) 


Where II is the photo generated current in Amperes 

Id is the ideal diode current in Amperes 

lo IS the leveise saturation current of the diode in Amperes 

q IS the electronic charge in Coulomb 

k IS the Boltzmann constant in Joule per Kelvin 

T IS the junction temperature in Kelvin 

Rg IS the series resistance in ohms 

Rah IS the shunt resistance in ohms 

n IS the ideality factor( the value of which depends on the recombination processes 
that are active m the cell) 
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Figure 4 3 1 V Characteristics of a solar cell 

An idealized curient voltage characteristic of solar cell is shown in figure 4 3 
One can note from the figure that the maximum output power of the solar cell 
IS lepiesented (A) the area of the maximum rectangle that can be drawn under 
the I V curve The point (A) usually leferred to as the "knee” of the I-V curve 
The electrical characteristics of solar cells are based on their I V curves The 
I V curve is based on the cell being under standard conditions of sunlight and 
cell temperature, and assumes there is no shading in the cell Since photovoltaic 
cells are electrical semiconductors, partial shading may cause the cell to heat 
up Undei this condition, the cells act as an inefficient conductor rather than 
an electiical generator Partial shading may rum shaded cells and also affect the 
power output of the cell Figure 4 4 shows the I-V characteristics of a shaded 
and unshaded cell 

The radiation and temperature are the mam factor which effects the maximum 
power point of a solar cell The effect of radiation and temperature on maximum 
power point of a solar cell is shown m figure 4 5 As indicated in figure 4 5, the 
maximum power point increases, if the radiation increases or if the temperature 
decreases [14] 


Some of the important definitions, which describes the cell properties found in 
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Figure 4 4 I V Characteristics of a shaded, and unshaded solar cell 


lOOOW/m 



Figure 4 5 Effect of temperature and radiation on I V Characteristics of a solar 
cell 
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literature [27] are presented here 


Fill Factor(FF) The ratio of the peak power to the product of the open circuit 
voltage ( Voc) and short circuit current ( Jsc) is defined as the fill factor (PF) 


FF = 


^raplmp 


^oclsc 


(4 4) 


where Imp and V-mp are the current and voltage at maximum power point 


Conversion Efficiency The conversion efficiency is the ratio of the optimal elec 
tiical power (Popt) delivered by the photovoltaic module to the solar mso 
lation (jBe)) received at a given cell temperature, T 


Popt 

~AEe 


(4 5) 


Where the optimal power Popt = VocIacFF is in Watts, Eg is in Watts 
per square meter and the cell area, A, is in square meters 


4 3 Forecasting Methodologies 

According to the basics of photovoltaic power generation and solar cell charac 
teristics, one can say that photovoltaic power generation mainly depends on the 
irradiation and the cell temperature It is very difficult to measure cell tempera 
ture, so instead of that, weather parameters like temperature and wind velocity 
are utilized In this study irradiation, temperature and wind velocity aie utilized 
foi forecasting photovoltaic power 

Here two forecasting methodologies are presented which were used by the author 
The first one is based on artificial neural networks The second one is based on 
multiple regression and statistical analysis 

4 3 1 Neural Network Based Prediction Model 

The artificial neural network used in this study, is a multi layered feed-forward 
network using a back propagation training algorithm The detailed explanation of 
this kind of network is explained in chapter 6 Here the configuration of proposed 
model for estimating the photovoltaic power and its training process is explained 
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Figure 4 6 Configuration of the neural network 


4 3 11 Configuration of Artificial Neural Network 

The configuration of proposed neural network is shown in figure 4 6 The network 
has three layers, namely input layer, hidden layer and output layer The input 
layer has four nodes for irradiation, temperature , wind velocity and a bias signal 
of 1 The hidden layer has eight nodes including one bias node The output layer 
has only one node, and the output signal is the predicted maximum power The 
nodes aie classified into two groups as white nodes and grey nodes In the white 
nodes no activation function is used, which means that the input goes directly to 
the output In the grey nodes unipolar continuous sigmoidal activation function 
(see chapter 6) is used The bias nodes are used here to improve the learning 
speed m the training process [12] 
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4 3 12 Initialization of The Training 


In all types of neural networks the training should be started by using proper 
weight matrices for output and hidden layeis Here the weight matrices are 
chosen randomly in between [0 5 0 5] except 0, for both output and hidden 
layers Before the training process the input vector should be normalized Here 
the input vector is normalized in between 0 1 to 0 9 by using the below equation 
where and Xnorm are the original and normalized input values 


0 8 


X 


norm 


max(X) — rmn(X) 


{xori — rmn{X)) +01 


(4 6 ) 


where X is consider the input vector This normalization is to give the 10 
percent variation of the minimum and the maximum range of the input vector 
The completion of the training process is done by checking the mean square error 
{MSE) This can computed by the equation below 


MSE = 


1 


Np I 


(4 7) 


where Np is the number of training patterns, I is the number of output neurons, 
in this case it is one di and yi are the desired and actual outputs The training 
process completes if the MSE will be less than or equal to the specified value e 
Here e is chosen as 0 001 If this is not met then one should give the maximum 
number of iterations {ITER), that the neural network should be trained The 
learning rate ( 77 ) and the momentum (/?) are chosen 0 5 and 0 3 During the 
training learning rate is changed as explained in chapter 6 , section 6 5 


4 3 2 Regression Based Prediction Model 

The regression model used in this study, is simple linear regression model with 
multiple inputs The detailed analysis of linear regression with multi variables 
is explained in chapter 5 The linear regression equation for photovoltaic power 
generation can be written as shown below 


(4 8 ) 


Power = a + 61 Rad + 62 Temp + 63 Ws 
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Where Rad is the irradiation Temp is the temperature and Ws is the wind 
speed 

By using regression analysis, the above regression coefficients a, 6i, 62 and 63 are 
to be found 

4 4 Data Analysis 

The data that is used for performing neural network and multiple regression 
forecasting models was taken form an experimental solar cell configuration which 
belongs to the Institut fuer Elektrische Energieversorgung, Technical University, 
Darmstadt, Germany In that configuration there are twelve modules and each 
module is having 36 solar cells 

The solai cells are of mono crystalline type and they are connected serially in 
this configuration These solar cells aie from Siemens and the specification type 
IS M55S Some main characteristics like peak power, short circuit current and 
open circuit voltage are presented below for one photovoltaic module 


Maximum Power (Pmax) 
Short Circuit Current {I sc) 
Open Circuit Voltage {Vqc) 


53 0 (W) 
3 27 (A) 
21 8 (V) 


Table 4 1 Specifications of M 55 S photovoltaic module 


In this configuration measurements weie taken foi each day with one second 
inteival The measurements are carried foi radiation, temperature and wind 
velocities, weie taken 24 hours for each day with one second interval The data 
used in this model are from Januaiy and June 1994 January data is for winter 
period power forecasting and June data is for summer period power forecasting 
In the summei periods, the data period is considered from 3 a m to 9 p m During 
the winter season the data period is considered from 8 a m to 6 p m The noise in 
the data was filtered and modified before it was used m the models Through the 
modification of the data, the values were reduced to one value for each minute 
This was done by taking the average for each minute interval data 
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4 4 1 Selection of Training Data Sets 

For the tiaming of the proposed neural network training the data sets were divided 
into two groups One is summer data set and othei is winter data set As 
explained in previous section, summer and winter data are having significant 
variations, so it is better to get two neural network structures separately for 
summer and winter data sets To get proper neural network structure m each 
season, three consequent days are selected m each season The training for each 
season is carried in two phases In the first phase the middle day is used to tram 
the net and other two days are used to test the network In second phase the 
first two days are used to tram the net and third day is used to test the net The 
data sets that aie used in summer and winter power forecasting are given below 

Summer data sets June 20th ,21st and 22nd of 1994, are selected to get proper 
neural net for summer period These figures are shown in Appendix A 5, A 6 
and A 7 In these figures power, radiation, temperature and wind velocity 
are plotted with respect to time In phase I 21st day data is selected to 
tram the net and 20th and 22nd days are used for testing In phase II 20th 
and 21st days are used to train the net and 22nd day data used for testing 

Winter data sets January 9th, 10th and 11th of 1994, are selected for winter 
power forecasting These figures are also shown in Appendix A 8, A 9 and 
A 10 In phase I 10th day data is chosen for training purpose, and 9th 
and 11th data’s aie used for testing In phase II combmedly 9th and 10th 
data is used to tiain the net and 11th day data set is used for testing the 
accuracy of the neuial net 

The training and testing results of both the neural networks for summer and 
winter are discussed in next section 
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Figure 4 7 Phase I 22nd june data neural network results 

4 5 Results and Comparisons 

4 5 1 Forecasting Results by Using Neural Networks 

Training of the proposed neural network model is earned for summer and winter 
seasons in both phases separately For the summer period for phase I the neural 
network results are shown in figuie 4 7 For the better view of the estimation 
accuracy two specific hours data was chosen for testing the proposed model, they 
are 14th hour of 20th june data and 9th hour of 22nd june data The results 
of these two hours aie presented in Appendix A 13 and A 14 Similarly for the 
second training phase results are shown in figure 4 8 and for better view, 9th 
hour of that day is chosen for testing and is shown in Appendix A 15 The 
results corresponding to the 20th june were also presented m Appendix A 12 
The accuracy test of the proposed model is done by calculating the absolute 
percentage error that is 


APE = 


\AC - ES 


X 100 


(4 9) 


AC 
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Figuie 4 8 Phase II 22iid june data neural network results 

Wheie AC is the actual measured value and ES is the estimated value of the 
powei 

In the similar manner it is carried for the winter data also The results for the 
winter season is shown in figure 4 9 for the phase I neural network For phase II 
the results are shown m figure 4 10 

4 5 2 Forecasting Results by Using Multiple Regression 

In Regression analysis the data is also divided into two seasons and for each 
season the analysis piocess is carried m two phases For summer period, in phase 
I, 21st June data is taken for regiession analysis and tested with other two days 
In phase II total two days that is 20th and 21st are chosen for regression analysis 
and IS tested for the third day The statistical analysis of the data is carried for 
20th June, 21st june and for both which are shown in Table 4 2 The value of , 
which IS called multiple regression coefficient, is nearer to 1, foi each day and the 
total of the both days, is very encouraging (see chapter 5) 

It means that the analysis is very good and it can estimate the other day power 
more accurately The forecasting results are shown in figure 4 11 for the phase I 
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Figure 4 9 Phase I 11th January data neural network results 



Figure 4 10 Phase II 11th January data neural network results 


I 


M 
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Figure 4 11 Phase I 22nd june data multiple legression results 



Figure 412 Phase II 22nd june data multiple regression results 
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Data File 

Intercept (a) 

bR 

bp 

bw 

B^ 

B 

20th June 
21st June 

0 0604 

0 05804 

0 8788 

0 82299 

0 11841 

0 04062 

0 05608 

7 198E-4 

0 946528 

0 94329 

0 97289 

0 971235 

Total 

0 05677 

0 8401 

0 08203 

0 05774 

0 94404 

0 971619 


Table 4 2 Statistical results for summer data 


and for phase II it is shown in figure 4 12 Again for the better view of the results 
hourly data is also chosen and is shown in Appendix figure A 18 and figuie A 19 
for phase I and figure A 20 for phase II 


Variable 

b 

Sb 

Fb 

Correlation with Power 

Radiation 

0 82299 

0 007645 

11589 95 

0 97079 

Temperature 

0 04062 

0 011439 

12 60846 

0 57665 

Wind Speed 

7 198Eh04 

0 02543 

8 013E-4 

0 033435 



Bpf^RaT) = 0 971235 

RaT) ~ ^ 000000 

Fr2 = 5966 74 


Table 4 3 21st June data multiple regression analysis 


The statistical results for 21st and total of 20th plus 21st are presented in Table 
A 1 and Table 4 4 The significance of B? and also the regression coefficients 
are tested by the F test (Fft^and Fb) and the results are presented in those ta 
bles, which are good Test of significance of B? and regression coefficients is 
discussed in chapter 5, section 5 8 The results of significance test are compa- 
rable with the null hypothesis value (see chapter 5) The standard error’s of 
the regression coefficients (Sb) and semi partial correlation coefficients (jR^p(TRa) 
andB^ P(WTRa) ) are also presented with Bp^tro) The standard error of re 
gression coefficients is less which means the regression equation is better fit for 
the historical data The semi partial correlation coefficients tells about the effect 
of each independent variables, that is radiation, temperature and wind speed, 
when they are presented in the regression analysis sequentially 
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Variable 

b 

Sb 

Fb 

Correlation with Power 

Radiation 

0 8401 

0 00597 

19790 8 

0 97006 

Temperature 

0 08203 

0 00787 

108 548 

0 62745 

Wind Speed 

0 05774 

0 015198 

14 4322 

0 360015 



Rp (RaT) = 0 971421 

^p(WiiaT) = 3 8472E-004 

Fr'^ = 11787 448 


Table 4 4 20th + 21st June data multiple regression analysis 



Figure 4 13 Phase I 11th January data multiple regression results 

Similar procedure is applied for the winter data also and forecasting results are 
shown in figuie 4 13 corresponding to phase I and figure 4 14 corresponding to 
phase II The results for the 9th January data were shown in Appendix figure 
A 16 From the statistical analysis the values of B? for all the days are very 
fax from 1 It means that the data is not good to fit the regression equation 
Although significant test results shows very fai from the null hypothesis value (see 
chapter 5) , the standaxd error of the legression equation not worth mentioned 
Statistical analysis of the winter data that is 9th, 10th and total (9th-t-10th) is 
also cairied which are shown in Table 4 5 Individual statistical results for 10th 
and total (9th+10th) are presented in Table 4 6 and Table 4 7 Statistical results 
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Figure 4 14 Phase II 11th January data multiple regression results 


for the 9th January is shown in Appendix A 


Data File 

Intercept (a) 


bx 

bw 


R 

9th January 
10th January 

0 00856 

0 02861 

0 91547 

0 25366 

0 12543 

0 30732 

0 00863 

0 03667 

0 6695 

0 42369 

0 81822 

0 65091 

Total 

0 01258 

0 61349 
1 



0 57422 

0 75778 


Table 4 5 Statistical Results for winter data 


Variable 

6 

Sb 

Fb 

Correlation with Power 

Radiation 

0 25365 

0 03212 

62 37185 

0 5118 

Temperature 

0 30732 

0 2377 

167 152 

0 60166 

Wind Speed 

0 03667 

0 02099 

3 0508 

0 06176 



Rp {RaT) = 0 64864 


= 146 0576 


Table 4 6 10th January data multiple regression analysis 
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Variable 

b 

St 

Ft 

Correlation with Power 

Radiation 

0 61349 

0 02795 

481 7477 

0 72705 

Temper atuie 

0 24722 

0 02235 

122 3128 

0 63455 

Wind Speed 

0 02646 

0 014185 

3 47998 

0 02444 


^%(TRa) ~ ^ 

RpiRaT) = 0 75695 

F-p(WRaT) ~ ^ 

Fr7 = 537 6658 


Table 4 7 9th + 10th January data multiple regression analysis 


4 5 3 Comparisons 

The compaiison study between neural network model and regression model is 
explained through the the following tables 

For summer period photovoltaic power forecasting, both the models gave better 
results in both phases The phase I and phase II comparison study results are 
shown in Table 4 8 and Table 4 9 The absolute percentage error for both models 
are below 10% which shows us the results are better and accuracy of estimation 
IS good 


Actual 

Powei(W) 

Estimated Power 
by regression 

Error 

(AC ES)/AC 

Estimated Power 
by neural net 

Error 

(AC - ES)/AC 

285 933 

273 1053 

0 04486271 

297 03376 

0 03882295 

248 11957 

242 82426 

0 02134174 

263 19864 

0 06077338 

313 97864 

304 26276 

0 0309444 

329 74286 

0 05020795 

367 62003 

342 24298 

0 06903064 

372 61959 

0 01359981 

421 95712 

392 04013 

0 07090055 

421 06262 

0 00211988 

281 7306 

275 64465 

0 02160198 

291 26263 

0 03383386 

521 8731 

516 5613 

0 01017839 

534 28768 

0 0237885 

339 97058 

327 42935 

0 03688916 

354 81731 

0 04367064 

33 43773 

30 891485 

0 07614882 

32 954305 

0 01445748 

33 177418 

31 382862 

0 05408967 

30 922437 

0 06796735 

31 31728 

31 255587 

0 00196994 

30 737674 

0 01850755 

31 113144 

31 867685 

0 02425153 

29 954991 

1 0 03722391 


Table 4 8 Summer forecasting results in phase I 
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Actual 

Power(W) 

Estimated Power 
by regression 

Error 

(AC ES)/AC 

Estimated Power 
by neural net 

Error 

(AC ES)/AC 

285 933 

268 78052 

0 05998781 

272 63153 

0 04651954 

248 11957 

241 03696 

0 02854515 

240 38665 

0 03116611 

313 97864 

298 01813 

0 05083311 

303 5089 

0 03334538 

367 62003 

339 5232 

0 076429 

352 07792 

0 04227764 

421 95712 

391 26404 

0 07273982 

401 64712 

0 04813286 

281 7306 

276 10715 

0 01996036 

268 8079 

0 04586901 

521 8731 

508 85815 

0 02493892 

500 78032 

0 04041745 

339 97058 

323 47888 

0 0485092 

332 11618 

0 02310317 

33 43773 

32 367237 

0 03201449 

32 195081 

0 03716307 

33 177418 

33 57368 

0 01194376 

32 115856 

0 03199652 

31 31728 

31 322279 

0 00015963 

30 441757 

0 02795654 

31 113144 

33 97892 

0 09210821 

31 88181 

0 02470552 


Table 4 9 Summer forecasting results in phase II 


For the winter period photovoltaic power forecasting, it is found from the results 
that the neural network model is better than the regression based model But 
comparing with the summer period forecasting results, the winter period fore 
casting results are found bad This is because of high noisy in the winter period 
data The results corresponding to winter data is shown in Table 4 10 and Table 






Actual 

Powei(W) 

Estimated Power 
by legiession 

22 45075 

38 281857 

45 07448 

63 420948 

39 700428 

69 09451 

56 287014 

74 07306 

49 0921 

64 68476 

35 268646 

55 965977 

46 01911 

71 24653 

53 43403 

73 73588 

102 58027 

96 9917 

123 82854 

105 58305 

178 26665 

117 01831 

171 75298 

121 07685 

69 930626 

73 8231 

58 056194 

76 871765 

67 82848 

78 269226 


Error 

(AC ES)/AC 

Estimated Power 
by neural net 

Error 

(AG - ES)/AC 

0 7051482 

20 607161 

0 08211705 

0 40702567 

41 997104 

0 06827314 

0 74039716 

42 920444 

0 08110784 

0 31598842 

52 437921 

0 06838333 

0 31762064 

45 463991 

0 07390412 

0 5868479 

32 934294 

0 06618773 

0 54819435 

47 656461 

0 03557981 

0 37994233 

51 044439 

0 0447204 

0 05447997 

94 89121 

0 07495652 

0 1473448 

113 72289 

0 08161005 

0 3435771 1 

165 59396 

0 07108841 

0 29505238 

155 61626 

0 0939531 

0 0556619 

72 018864 

0 02986157 

0 3240924 

61 180204 

0 05381011 

0 15392858 

69 057034 

0 01811265 


Table 4 10 Winter forecasting results in phase I 


Actual 
Power ( W) 

Estimated Power 
by regression 

Error 

(AC ES)/AC 

Estimated Power 
by neural net 

Error 

(AC ES)/AC 

22 45075 

25 763172 

-0 1475417 

23 754934 

0 05809088 

45 07448 

60 269527 

0 33710983 

45 785301 

0 01576992 

39 700428 

61 46312 

-0 5481727 

42 10977 

0 06068805 

56 287014 

73 18864 

-0 3002757 

53 665741 

0 04656976 

49 0921 

62 826622 

0 27977055 

50 609778 

0 03091491 

35 268646 

46 833687 

0 3279128 

35 61555 

0 00983603 

46 01911 

65 174706 

0 416253 

46 731944 

-0 01548996 

53 43403 

70 86716 

0 32625517 

51 055351 

0 04451618 

102 58027 

117 29827 

0 1434779 

103 03717 

-0 00445408 

123 82854 

132 405 

0 06926078 

126 02267 

0 01771912 

178 26665 

151 21524 

0 15174687 

181 4929 

-0 01809787 

171 75298 

162 0112 

0 05671969 

180 92085 

0 05337823 

69 930626 

74 128975 

0 06003591 

65 69936 

0 06050662 

58 056194 

73 20502 

0 2609338 

57 035314 

0 01758434 

67 82848 

76 70768 

0 1309066 

63 767411 

0 05987263 


Table 4 11 Winter forecasting results in phase 11 
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4 6 Conclusions 

In this chapter photovoltaic power forecasting models weie discussed Two mod 
els weie constructed for photovoltaic power forecasting, a neural network based 
model and a multiple regression based model In neural network model standard 
multilayer feed forward back propagation learning algorithm was used Photo 
voltaic power forecasting is carried for two seasons, summer and winter Results 
were presented for both the models corresponding to summer and winter From 
the results it is concluded that the summer period power forecasting gave better 
results than the winter period power forecasting It is found due to the high noisy 
in winter data comparing with summer data For the summer period both the 
models gave good results and for the winter period neural network is better than 
the regression model The used data included recognizable noise and it would be 
interesting to test the models with more accurate data 


5 Multiple Regression 


5 1 Introduction 

The problem of finding the most suitable form of equation to predict one variable 
from the values of one or more other variables is called the problem of regression 
[7] In 01 del to estimate the regression curve of y on a: we must first specify the 
functional form of the cuive Some examples are 


Y = a + bX 

(5 1) 

Y = 0 + biX + b2X^ 

(5 2) 

Y = ae^ 

(5 3) 


Here only linear regression is considered with one Y variable and several X 
variables The analysis of linear regression model for multiple X variables is 
known as Multiple Regression Analysis Multiple Regression equation can be 
written as shown in the below equation 

Y' = a + biXi + b2-^2"f"» 

One h we to find the values of a, 6 i , 62 , M that will result in the highest 
possible positive correlation between the observed values Y and Y' This can be 
done by using historical data 
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5 2 Multiple Regression Analysis 

With k different X variables the multiple regression equation is 

y'' Q- -{• bxXx + 62 -X^ 2 +) i-k-b^Xk (5 5) 

The total sum of squared deviations of the Y values from the mean of the Y 
values {S Slot = can be partitioned into two orthogonal components 

the sum of squares for the linear regression (SSreg = ~ ^lie 

residual sum of squares (SSrea = — Y'Y) [8] 

SStot — SS^eg “b SSyes (5 6) 

= + (5 7) 

where Y is the mean of Y values To find the values of a, 61,62, ,6^, the 

residual sum of squares should me mininmzed Then the required value of a will 
be 


a = Y- 61X1 - - biXk (5 8) 

Where X is the mean of X values Substituting this value of a m 5 1, then 

y' = y + 61 (Xi ~ Xi) + 62(^2 - X2) + + b^{Xk - Xk) (5 9) 

Then, if X](y — X is to be minimized, the values of 6 must satisfy the following 
equations 

6 i + 62 ^(xixa) + +bk ^(xi®fc) = X^(®i2/) 

61 ^(X2®l) + 62 ^(®2a3fc) = J2(®2y) 

(5 10) 


61 X](a3fc®i) + b, J2(xkX2) + + bfc = Y^i^kV) 
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wheie 

Xi = X-i X ^5 11 ) 

y, — Y^-Y (5 12) 

The values of b can be obtained by solving the above equations After getting the 
values of b they can be substituted in 5 8 to find a The meaning of a regression 
coefficient in a multiple regression equation is If all variables except Xt are held 
constant, then 6t is the amount by which Yi increases with unit increase in Xi 
[ 8 ] 

5 3 The Correlation Coefficient 


The correlation coefficient between two variables Y and X can be defined as 




(5 13) 


The significance and the importance of the coirelation coefficient can be ex 
plained by analyzing the equation 5 2 

In case of two variables Y and X, if there is a perfect linear relationship be 
tween Y and X so that all of the plotted points fall precisely on a straight line, 
then SSres will be equal to zero and SSreg will be equal to SStot If there is 
no tendency for the Y values to be linearly related to the X values, then the 
legiession coefficient B and SSreg will be equal to zero and SS,e8 will be equal 
to SStot An important index of the degree to which the Y and X values are 
lineally related is obtained by dividing both sides of the above equation by SStot 
[ 8 ] 


putting 


then 1 

SSreg 

SStot 

then 1 


SS 


reg 


SStot 

J2 


+ 


SSr 


S S tot 

B^y, 4 " (1 ByYl) 


(5 14) 

{5 15) 
(5 16) 


where Byyi is the corielation coefficient oi in case of multiple regression it is 
called Multiple regression coefficient Obviously B^y cannot be greater than 1 
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Figure 5 1 The significance of Correlation Coefficient 
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and consequently, Ryv can range only within the limits 1 to 1 The significance 
of Ryy will be explained in case of one input (^) and one output (y) below 
[ 8 ] 

In case of two variables that is X and Y the significance of r is explained by using 
four giaphs shown in figures a,b,c,d Figure 5 la and figure 5 lb are the plots 
where r = 1 00 and r — 1 00, respectvly In figure 5 Ic, r value approximately 
0 8 and in figure 5 Id it is close to 0 3 The value of provides an index of 
the degree to which a set of plotted points cluster about the regression line The 
closer the points fall along the regression line, the larger the value of and the 
gieatei the piopoition of the total sum of squares accounted for by the linear 
regiession of Y on X When the value of is small as shown in Figure 5 Id, 
the plotted points will show considerable scatter about the regression line, and 
the proportion of the total sum of squaies accounted for by the linear regression 
of Y on X will be small 


5 4 Regression Analysis with Standardized 
Variables 


If all the variables in multiple regression problem that is k vaiiables of X and 
Y are in standardized form, then the value of a in the regression equation will 
be equal to zeio [8] In this instance all of the regression coefficients will be in 
standardized form, and the regression equation will take the form 

Zy' ~ bxZi + + ”b (5 17) 


The definition of standardized variable is one that has a mean value equal to 
zero and a variance and standard deviation equal to 1 Any variable Y can 
be transformed into a standardized variable by means of the following linear 
transformation 


;£ = 


Y -Y 
sy 


(5 18) 


where sy is the standard deviation of Y values 

If the residual sum of squares is to be minimized, then the regression coefficients 
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must satisfy the following equations 

biTii + 62^12 + + ~ 'i'Yl 

biVzi + b2r22 + + bkr2K = /’r2 

(6 19) 

bk'f'kl + b2rk2 + + bkTkk ~ Tyk 

wheie Tij IS the inter correlation coefficient between and £Cj, and Tyi is the 
coirelation coefficient between Y and x^ for 1 = 1,2,3, k The k X k inter 
coiielations of the X vaiiables can be represented by a symmetric matrix which 
is designate by Rtj The k correlations of the X variables with the Y variable 
can be represented hy & k X 1 column vector that is designate by Ryi Similarly, 
the k values of the regression coefficients can be represe noted by a ^ X 1 column 
vector that is designate by To find B, find the inverse of Rij and multiply 
with Ryi 

R-^RY^ = B, (5 20) 

One can obtain the original multiple regression coefficients that is B by using 
the below formula 

b. = i-6. (5 21) 

Sx 

From these values of b, one can obtain a by using 5 8 
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5 5 Mean Squares 

When a sum of squares is divided by its degrees of freedom (d f), the resulting 
value IS called a mean square The number of degrees of freedom can be thought 
of as the number of independent comparisons available [7] From the literature [8] 
the degrees of freedom of each sum of squares that is SStot, SSreg, and SSres, 
IS n 1 , k and n k 1 Where n is the number of samples, k is the number of X 
variables Then the mean squares of each sum of squares is 

MStot = (5 22) 

n — 1 

MS,,, = (5 23) 

MS ,, = - (5 24) 

n ~ k — 1 

5 6 Standard Error of the Regression Coefficients 


The standard error of a regression coefficient 6t when several X variables are 
involved will be given by 


Sb 


MSr 




(5 25) 


where is the squared multiple correlation of X, with the remaining X van 
ables When the inverse of the correlation matrix between the X variables is 
available, the value of is easily obtained It will be given by 



where r** is the diagonal entry in the inverse of the correlation matrix 


(5 26) 


5 7 The Semi Partial Correlation Coefficient 

A semi partial correlation coefficient is the correlation coefficient between two 
variables X^ and Xj after the variance that Xfc has in common with X. and 
Xj IS removed from only one of the two variables In case of three variables, the 
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semi partial correlation coefficient between Xi and X2 with the variance of Xs 
removed from X2 will be given by 


’’1(2 3) 


’’12 — ’’13’’23 




(5 27 ) 


23 


where the notation 1*1(2 3) indicates that it is the corielation coefficient between 
Xi and X2 aftei the vaiiance that X3 has in common with X2 has been elim 
mated from X2 Applying this for one Y and two X variables, semi partial 
corielation coefficient can be written as 


’’y-(i 2) 


ryi — rY2f’i2 


^A■ 


(5 28 ) 


' 12 


The relation between semi partial correlation and correlation coefficient is 


„2 _ D2 _ 

^y(i 2) ~ •”'yi2 'y2 

Genei aiming this for any number of X variables, 
be written as 


(5 29 ) 


then correlation coefficient can 


•^^123 A- ~ '^Vl ^y(2 1) +’’y(fcl23 jfc-l) (5 30) 

The theoretical behavior of above equation can be explained for three X variables 
as follows The proportion of the total sum of squares accounted for by Xi, 
given that it is entered first in the regression equation, will be After Xi 

IS entered in the regression equation, the increment m the proportion of the 
total sum of squares accounted for by X2, given that it is entered second in 
the regression equation, will be ry(2 1) Now that both Xi and X2 are in the 
regression equation, the increment due to X3 will be given by 1^1 


5 8 Test of Significance 

One of the most important thing m the multiple regression analysis is "test of sig 
nificancc” of multiple regression equation coefficients and correlation coefficient 
Significance tests for these values are carried by using statistical tests Commonly 
known tests m the statistics are F test, t test and test [7] 
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General procedure for making some test is explained here with an example [7] 
Assume theie aie n samples in the sample space which are taken from a normal 
distribution with unknown mean ^ The sample mean is denoted by x If one 
IS interested in a particular value and asks some question ’does x differ sig 
nificantly from then test of significance can be earned by calculating the 

observed ”t test’ value It can be found by using below equation 


to = 


X — [Xo 
sf\/n 


(5 31) 


Where s stands for the standard deviation of the samples Now the percentage 
point ta /2 o IS chosen so that there is a probability a of getting a larger observa 
tion from a t distribution with v degrees of freedom If |to| is larger than to- 
then the result is significant at the lOOo: percent significance level The values 
of tot V IS available in most of the statistical books [7] 

Similar test pioceduie is taken foi the correlation coefficient and regression co 
efficients As the above procedure is explained for the t test, same procedure is 
also carried for the F test also But the use of F test is it compares more than 
one sample space In this case the chosen percentage point can be denoted by 
■Fa /2 ui v 2 , where vl and v2 are the degrees of freedom for the two samples 


5 8 1 Test of Significance of Ryy' 


It IS often useful to perform a significance test to see if the observed Ryy is 
significantly difi!eient from zeio In other words it is desired to test the null 
hypothesis Ho Ryyi = 0 against the alternative hypothesis Hi Ryy/ ^ 0 
If all the regression coefficients B are equal to zero, then SSyeg will be equal to 
zeio and Ryy will also be equal to zero The obtained values of h[s are estimates 
of the corresponding population values Pt A test of the null hypothesis that the 
all the population values are equal to zero will be given by 


Fo 


R^y /fe 

(1 ~ R^Y,)/(n - k-1) 


(5 32) 



[ 
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wheie k the number of X variables If both the numerator and the denominator 
of eq aie multiplied by SStoti then the equivalent test, 

F SS,eo/k ^ MS, eg 

° SS,es/{n-k-l) MSrs 
For both F latios the degrees of freedom for the numerator will be equal to k, 
the number of X variables, and the degrees of fieedom for the denominator will 
be equal to n ^ i where n is the number of samples This observed value is 
compared with the Fa /2 m v 2 , with different a values Here the values of vl and 
V are k and n k 1 


(5 33) 


5 8 2 Test of Significance of B 

For a test of significance of a regression coefficient in a multiple regression one 
cm find the observed value of the test by using below equation 

= , *’• (6 34) 

The square of t test will be equal to F test with 1 and n ^ i degrees of freedom, 
so one can do F test for this and compare with the Fa /2 vi v 2 Iii ^tiis case the 
values of and v are 1 and n k 1 


5 9 Multicollinearity 


The existence of high inter correlations among the X variables is called the prob 
lem of mull icolhnearity Evidence regarding a high degree of multicolhnearity 
c m be provided by the standard errors of the regression coefficients Form below 

equation 


« - (5 35) 

wheie IS Ihe squared multiple correlation of X, with the remaining X van 
ables As approaches 1, the denominator becomes small and Ss becomes 
laige An extremely large value of Ss for any of the regression coefficients may 
be an indication of multicolhnearity 
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5 10 Algorithm for Multiple Regression 


The detailed steps that have to be done m Multiple Regression Analysis are given 
IS this section Multiple regiession equation, as explained before is given below 

Y' = a + biXi + 62-^2"h> » 4'hfcX^fc (5 36) 


One have to find the regression coefficients 0,61,625 j h* by using the his 
toiical measuiements between independent variables (X) and dependent variable 
(y) Aftei that one have to do statistical analysis of the data and test of signif 
icance of the 1 egression coefficients The procedure of all analysis is carried from 
Step 1 to Step 12 


Step 1 Load the hisioncal data of Independent (Xkxn) depen 
dent vanables (1 X N), where N is the number of points in 
the tiaming data bank and K is the number of Independent 
vanables 

Step 2 Form the Intercorrelation (Rij) and Coirelation matrices (R^J ) 
Step 3 Find the inverse of Intercorrelation matrix that is R~/ 

Step 4 Find B, which is called standard form of regression coeScient 
vector by below equation 

JB, = 

Step 5 Now £nd the regiession coefficient vectoi (B — [bi^bz, bk]) 

by using below equation 

b. = ^b, (5 38 ) 

SX 

wheie i = 1,2,3, ,k , sy and sx are the standard 
deviations of dependent and independent vanables 
Step 6 Calculate intercept coefficient (a) by using below equation 

a = y - 61 X 1 , , -^kX k (5 39) 


where Xt and Y are the mean values of independent and 
dependent variables 
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Step 7 


Step 8 


Step 9 


Step 10 


Calculate the multiple correlation coefRcient ) by below 

equation 

B^y = (5 40) 


Now find all sum of squares which are denoted by SStot , 
SS,eg and SS,ea by using below equations 

SStot = 

(5 41) 

wheie y — Y — Y 


SS,eg = ^ + ^2 a! 22 / + 

+ it 52 

(5 42) 

where x, ~ X, — X 


SSrea — SStot SSreg 

(5 43) 


Find the Mean Squares of each sum of squares by using below 
equations 


MStot 
MS,,, 
IVf S, ea 


SStot 
Sfofi], f 
SSreg 

Sregdf 

SSrea 

Sresdf 


(5 44) 
(5 45) 
(5 46) 


wiiere d f stands for degrees of fieedom 

The degiees of fieedom for each sum of square is N - 1 , K 

and N — K — 1 

Find the siandaid error of all regression coefficients that is 


B' s by below equation 


M5, 


rea 




(5 47) 


wheie = 1 — m which r” is the ith diagonal element 
in B~ 


1 


6 Artificial Neural Networks 


6 1 Introduction 

Aitifici'il neural networks are parallel models comprrsed of densely rnterconn 
ectcd rdaptive processnrg unrts These networks are line grained parallel imple 
mentations of nonlinear static or dynamic systems A very important feature of 
these networks is then adaptive nature, where ’’learning by example” replaces 
’’programming” in solving the problems This feature makes such computatronal 
models very appealing m application domains where one has little or incomplete 
understanding of the problem to be solved but where training data is readily 
available [10] 

Artificial neural networks are viable computational models for a wide variety of 
problems These include pattern classification, forecasting and prediction, clus 
teiing, nonlinear system modeling etc In this chaptei it is mainly concerned on 
the ’’forecasting and prediction” type application by using the well known Error 
Back Pi opagaLion (EBP) algorithm 

1 his chaptei is organized mainly in five parts in winch, section 6 2 discusses 
the basic model of the artificial neural networks In section 6 3 consists about 
the learning process in the neural network concentrating mainly on Gradient 
Descent Based learning rule In section 6 4 detailed description of EBP algo 
iithm IS explained, where as section 6 5 points the modifications of EBP and its 
limitations 
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Figiue 6 1 The basic model of a neuron 


6 2 The Basic Model of Artificial Neural Network 


The fund xmcntal to any neuial network is it’s information processing unit, which 
IS called is new on The block diagram of figuie 6 1 shows the model of a neuron, 
which foims the basis for designing artificial neural networks [11] Every neuron 
consists thiee basic elements which are called synaptic connections, summing 
point and an activation function, which aie shown in figuie 6 2 

Synaptic connections Synaptic connections are mainly the connections between 
iiipui signal (»t) of tli-O neuion to the summing point It is characterized 
by 1 weight {w) or strength of its own As an example a connection weight, 
between ilh input signal (a;,), to thejth summing point can be denoted by 

Summing point It is the point wheie all the input signals multiplied by their 
synaptic weights will be added 

Activation function This is for limiting the amplitude of the output [yk) of a 
neuion The activation function is also lefeired to as a squashing function 
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in that it squashes (limits) the permissible amplitude range of the output 
signal to some finite value Typically, the normalized amplitude range of the 
output of a neuion is wiitten as the closed unit interval [0 1] oi alternatively 

[l,i] 

The neuional model of figure also includes an externally applied connection which 
IS called a bins This is denoted in figure as bj~ The bias has the effect of 
incieasing oi loweimg the net input of the activation function, depending on 
whethei it is positive or negative, respectively In mathematical terms, neuron k 
can be desciibed by wiitmg the following equations 

m 

= + (6 1) 

1=1 

and 

Vk = finetk) (6 2) 

wheie a-ie fhe input signals, Wfci,‘uifc 2 , Wkm are the synaptic 

weights of ncuion k, netk is the linear combiner output due to the input signals 
and bias, 6^, is the bias signal, /( ) is the activation function, and yk is the 
output signal of the neuion 

6 2 1 Types of Activation Function 

Mostly in neuion models choosing the activation function is important Here 
various activation functions which aie found in literature [11] will be explained 


Threshold Function The graphical form of this function is shown in figure 6 2 
(a) Ihe mathematical lepiesentation of this function is 


finet) 


1 if net > 0 
0 if net < 0 


(6 3) 


Such a neuion is referred to in the literature as the McCulloch Pitts model 
In this model, the output of a neuron takes on the value of 1 if the net of 
the neuion is nonnegative, and 0 otherwise 
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Piecewise Linear Function Tbs is shown m figure 6 2 (b) , the mathematical 
lepiesent ation of this function is 

{ 1 if net > +1 /2 

net if net < +1/2 > net > -1/2 (6 4) 

0 if net < —1/2 

whcic the amplification factor inside the linear legion of operation is as 
sumecl to be unity This form of an activation function may be viewed as 
an appi oxtmation to a nonlinear amphfiei 


Sigmoid Function The sigmoid function, whose graph is s-shaped, is by far 
the most common form of activation function used in the construction of 
aitifici il ncuial networks It is defined as a strictly increasing function that 
exhibit b a balance between lineai and nonlinear behavior An example of 
tlie sigmoidal function is the logistic function, difined by 


finet) 


1 

1 + exp{—\net) 


(6 5) 


whole A > 0 is called the steepness parameter of the sigmoidal function 
It’s giaphical view is shown m figure 6 2 (c) for different values of A The 
main feature of this type of activation function is, it is differentiable, where 
as threshold function is not In Error Back Propagation which is discussed 
latei in this chapter, is used maily sigmoidal function, which uses its differ 
entiable 1 unction 


riic iDias IS an external paiametei of aitificial neuion k, can be formulted as 
a anotlui input signal aio, which has a value of il, connected with a synaptic 
weight of value b Now the model of a k th neuion can be written by following 
equations 

m 

— and 

Vk = fineth) 

wheic X IS input signal vector, having xq = il, "w is the weight vector having 

Wfeo = h 
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As explained before that neural networks are densily interconnected neurons it 
IS possible to get diffeicni types of neural network stiuctuie by choosing different 
neuion models with diffeient activation functions One of such type of network 
lb Backpiopagaiwn neural network which is explained later The other mam 
impoitant thing m iieuial network is ’’learning algorithm” It is explained m 
comming section 

6 3 Learning Rule 

One of the most significant attributes of a neural network is its ability to learn 
by interacting with its environment or with an information source Learning in a 
neural network is noiin<illy accomplished through an adaptive procedure, known 
as a leaining luk or algorithm, whereby the weights of the network are mere 
mentally adjusted so as to improve a predefined performance measure over time 
The learning process can be viewed as ’’search” in a multidimensional parameter 
(weight) space for a solution, which gradually optimizes a prespecified objective 
(ciiteiion) function 

There arc number of learning rules presented in the literature [10] which are 
basically classified into three types of learning tasks They are supervised, un 
supervised and reinforced learning tasks In supervised learning (also known as 
learning with a teacher), each input pattern leceived from the environment is 
associated with a specific desired target pattern Usually, the weights are syn 
thesizcd gradually, and at each step of the learning process they are undated so 
that the Cl 101 between the network’s output and a corresponding desired target 
is reduced On the other hand, unsupervrsed learning (known as learning with 
out teacher) involves the clustering of unlabeled patterns of a given training set 
Reinforcement learning involves updating the network s weights in response to 
an ’’evaluative” tcachci signal, this differs from supervised learning where the 
teachei signal is the co? i ect answei 

The following genei al learning i ule is adopted in neural network studies [28] The 
weight vector = [Wti , Wiz ? j '>^tnY increases in proportion to the product 
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Figuie 6 3 Illustiation foi weight learning rules 

of input \ and leaining signal i The learning signal r is in general a function 
of wit, ai, and sometimes of the teachei’s signal d, The whole thing is shown in 
hguie 6 3 

r — r{wt,x,dt) (6 7) 

The increment of weight vectoi wi, produced by the learning step at time t ac 
cording to the general learning rule is 

Awit(i) = cr[wi(t),x{t),d^{t)]x{t) (6 8) 

where c is a positive number called the learning constant that determines the rate 
of learning f he weight vectoi adapted at time t becomes at the next instant, or 
learning step, 

w^{t + 1) = Wi{t) + cT[w,{t)^x(t)^dt{t)]x{t) (6 9) 

T-his can be written in formal notations 

^ cr{w^,x’^ 4)®*= (6 10) 

where is the weight vector at k th step and is the k-fl th iteration 

weight vectoi 
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Figure 6 4 Delta Learning Rule 

In [28], vaiioufa learning algorithms aie discussed, based on this basic learning 
algoiithm Among them ’’Gradient Descent based Delta Algorithm” is one which 
gives basic idea about the Error Back Propagation Algorithm The mam aim of 
this chaptei is to explain the EBP, so it is better to know the ’’Delta Rule”, 
whidi IS explained latei m this section 

6 3 1 Delta-Rule 

The delta learning lule is only valid foi continuous activation functions like figure 
6 2 (c) lh( Icainmg signal foi this is called delta and is defined as [28] 

7 ^ [dt -- f {net)] f {net) (6 11) 

lh( tcim f'{net) is the deiivative of the activation function f{net) computed 
foi net = wlx Ihc pioceedmg of the delta learning lule is indicated in figure 
6 4 1 Ills k ammg rule can be readily derived fiom the condition of least squarred 
euoi between j/i and d, The sum of squaired erior defined as 
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where i/t = f(wlx) The gradient of this error is given by 

VE = -{d,~y,)f{w\x)x 

The conrponcnlb of gr idieni vector arc 

dE 

~{d^~y^)f{w\x)x for j=l,2, ,n 


dw 




(6 13) 


(6 14) 


Since ihc niinimiJ’ation of the error requires the weight changes (AW) to be 
the negative gradient direction, so 


in 


Aw, = -r)VE (6 15 ) 

where rj is a positive constant known as learning rate Then one can obtain from 
equation 6 1 ] and equ li ion 6 15 

Axv,j = Ti{d, ~ y,)f{net,)xj for j = l,2, ,n (6 16) 

One disadvaiitexge of the delta learning rule is , it progress very slowly when the 
f’(nct)=0, even when the eiroi is large f’(net) = 0 when net has large magnitude 
To avoid this the common technique is replace /' by f plus a small positive bias 
e [10] In this case the weight adjustment becomes 

Aw,j — ri{d, - yt){f{net,) + e)xj for j = 1, 2, , n (6 17) 

one of the piimaiy advantages of the delta rule is that it has a natural extension 
that may lx use d to tram multilayered neuial nets 


6 4 Back Propagation 

Ihe back propagation algorithm is a supervised learning method for multi layered 
feedforward neural networks It is essentially a gradient descent local optimization 
technique, it involves backward error correction of the network weights The 
backpropagation learning rule is central to much current work on learning in 
artificial neural netwoiks Back propagation trained multilayer neural nets have 
been applied successfully to solve some difficult and diverse problems, such as 
time scries prediction, non linear system modelling, pattern classification etc [10] 



Chaptei 6 Aitificial Neuial Netwoiks 


Page 71 


-< dj 

y, d, 

Yl ^ 

Figuie 6 5 A two layci fully intei connected feedforward neural network archi- 
ItciuK 

6 4 1 Learning Rule for Multilayer Feedforward Neural 
Networks 

Consider the two 1 lyci feedforward architecture shown in figure 6 5 This network 
receives a set of sc ilai signals Xo, aij, , Xn where ico is a bras signal equal to 
i This set of signals constitutes an input vector x € The layer receiving 

the input signal is called the hidden layer Figure 6 5 shows a hidden layer having 
J units 1 he output of the hidden hyei is a (J+l) dimensional real valued vector 
z — [zo^Zxi ^zj] Again, Zo = -1 represents a bias input The vector z 
supplies the nipul for the output layer of L units The output layer generates 
an L dimensional vcctoi y iJi i espouse to the input x which, when the network 
IS fully 1 rained, should be idcuiical (or very close) to a desired output vector d 
associated wilJi x 

The activationfunc tion fh of the hidden units is assumed to be a differentiable 
nonlinear function (typically, fn is the logistic function defined by fh{net) = 
1/(1 with values of A close to unity) Each unit of the output layer 

IS assumed to have the same activation function fo The functional form of fo 
IS determined by the desired output signal/pattern representatron or the type of 
application I’or example, if the application is of estimation type with desired 
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values leal positive, then a satuialiiig nonlmeaiity similar to fh may be used for 
/, rurally denote w,, the weight of the j th hidden unit associated with the 
input signd Siinihily vi, is the weight of the 1 th output unit associated 
with the hidden signal Zj 

Next, eonsidc i x set of m input/output pairs where df. is an L dimensional 

vector representing the desired network output upon presentation of Xk The 
object here is to xdaptivcly xdjust the w, v weights of this network such that the 
underlying function/mapprng lepiesented by the training set is approximated or 
learned Since the learning here is supervised, and error function may be difined 
to measure the degree of approvimation for any given setting of the network’s 
weights A conimonly used eiioi function is the sum of squaired error (55 E) 
measure Oiuc a suitable enoi is formulated, learning can be viewed as an 
optimi/ation pioceas Ihat is, the eiror function serves as a criterion function, 
and the Icaiiiing algorithm seeks to minimize the criterion function over the space 
of possible weight settuxgs If a differ entiable criterron function is used, agradient 
descent on such a function will naturally lead to a learning rule Now consider 
the following eiioi function which is minimized by adjusting the weights 

(6 IS) 

Here w represents the sot of all weights in the network 

6 4 2 Error Backpropagation Learning Rule 

Smcc the largois for the output units are explicitly specified, one can use the 
delta lulc dncctly for updating the vij weights That is 

&E 

■where Si = (di — yi)foineti) (6 20 ) 

for L and j = 0,1,3 J 

Here neti = VijZj is the weighted sum of the 1 th output unit, /' is the 

dciivative of fo with respect to net, Si is the error signal of the 1 th output 


Eiw) = - yiY 


1=1 
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unit and and lepieseids the updated (new) and current weight values, 
icspcc lively T-hc Zj values aie computed by propagating the input vector x 
tluougli the hidden liyei according to 






fh^nstj) j — 152 , J 


(6 21 ) 


a=0 


dhe le lining rule for the hidden layer weight Wjt is not as obvious as that for 
the output layer because a set of target values for hidden units are not avail 
able However, one may derive a learning rule for hidden units by attempting to 
minimise the output layer error This amounts to propagating the output errors 
back thiough the output layer toward the hidden units in an attempt to estimate 
dynainu taigrts for these units Such a learning rule is termed error backpropa 
gation or tin hadpropagation learning rule and maybe viewed as an extension 
of tin d( It r inh 

lo comph t( the derivation of backpiop for the hidden layer weights, and simi 
hr to the pi( coding derivation for the output layer weights, gradient descent is 
peifoimc cl on the criterion function in equation 6 18 This time the gradient is 
calculate d with respect to the hidden weights 

Awn = — T? --— where j = l,2,3, = 0,1,2, n (6 22) 

dwji 

Using the chain rule for differentiation, one may express the partial derivative in 
equation 6 22 as 

dE dE dzj dnetj (6 23 ) 

dwj, dzj dnetj dWjt 



with 


dnct, 

div,, 

9zj 

dnetj 

ind 

dL 

dzj 


'^ I'j 


d 1 ^ 

~ fo{neti)f} 

^ 1=1 

- EW ~f- 

t=0 OZj 

L 

t=0 

L 

- 8iwij 
t=0 


combining c qtial ion 6 25 uid equation 6 26 then 


(6 24) 
(6 25) 


(6 26) 


L 

(6 27) 

1=0 

whcie Sy IS tlu c 11 01 signal of j th hidden neuron The weight adjustnaent m the 
hidden layci becomes 

Aiiiji = TfSyXt (6 28) 

It IS usually jxissibh to oxpiess the deiivatives of the activation functions in 
equation 6 19 and (qiutiou 6 27 in teims of the activations themselves For 
cxamjik, foi tin logistic activation function, f{net) = 1/(1 

f{nei) = Xf{net)[l — /(net)] (6 29) 

and foi the liypoibolic tangent activation function f(net) = tanh(l3net) 

f{nGt) = /3[1 •- /^(net)] (6 30) 

These Icainmg equations may also be extended to feedforwaid nets with more 
than one hidden layei and/oi nets with connections that jump ovei one or more 
layei s 
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6 4 3 Error Back-Propagation Training 

rigUK 6 6 illusiiatcs ihc flowchait of the eiroi back propagation training algo 
iitlun loi 1 1) iMC two 1 lyu network as m figure 6 5 The learning begins with the 
fee dloi waul Kcall phase (Step 2) After a single pattern vector x is submitted 
it ihc uipul ) the layers responses x and y are computed in this phase Then, 
the error srguil coixiput ition phase (Step 4) follows Note that the error srgnal 
veetor riiiisl be determined in the output layer first, and then it is propagated 
tow lid the network input nodes The L X J weights are subsequently adjusted 
within lire matrix V m Step 5 Pinally J X N weights are adjusted within the 
m ill IX W in Step 6 Note that the cumulative cycle of input to output mapping 
IS computed 111 Step 3 as r sum over all continuous output errors m the entrre 
training sc 1 Ihc final uroi value for the entire training cycle is calculated after 
each ( oinplc led pass thiough the training set Si, 012 , Xp The learning proce 
duK stops when the final error value below the upper bound, Emax, is obtained 
as shown in Step 8 [28] 

Figuu 6 7 depicts lire block diagram of the error back propagation trained net 
work operation and shows both the flow of signal, and the flow of error within 
the network [28] lire feedforward phase is self explanator The shaded portion 
of the diagram 1 efers to the feedforward recall The blank portion of the diagram 
lefei s to the training mode of the network The back propagation of error d-y 
from each output, for f = 1, 2, 3, L, using the negative gradient descent tech 
mque IS divided into functional steps such as calculation of the error signal vector 
Si and ( alculation of the weight matrix adjustment AV of the output layer The 
diagram rho illustiates the calculation of internal error signal vector Sj and of 
the resulting weight adjustment AW of the input layer 

The comjrlete procedure for updating the weights in a feedforward neural net is 
summui/cd below for the two layer architecture of figure 6 5 
Given there are P training patterns 

Step 1 Set the learning latc y to a small positive value Initialize 
all weights at small random values and refer them as current 
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rigiuc 6 6 Eifor Back Propagation tiammg {EBPT) algorithm flowchart 
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Figuie 6 7 Luoi Back Piopagation tiainmg (EBPT) algorithm block diagram 


weights Choose Emax, p = 1,E = 0 
Step 2 Select an input pattern Xk from the training set (preferably at 
1 indom) and propagate it through the network, thus gener- 
ating hidden and output unit activities based on the current 
weight settings 

Step 3 ('omputc the eiioi E ~ E -{■ l{dk — VkY 
Step 4. Compute the erioi signals Si and <5, by employing the equa- 
tions 6 20 and 6 27 

Step 5 Adpxst output layer weights V as shown in equation 

^new ^ ^oid (6 31) 

Adjust hidden layei weights W as shown m equation 

^new „ ^ (6 32) 


step 6 
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Step 7 if p < P then p — p 1, md goto step 2 other wise goto 
step 8 

Step 8 Oiu ti unwg cycle is completed Now test for convergence 
]l E <. Emax ihen teimmate the training session, else put 
P = 0 ind p ~ 1 initialize the new training cycle by going 
to step J 

It bhould be posbiblc that backpiopagation may fail to find a solution that passes 
the coiiveigence test In this case, one may try to reinitialize the search process 
tune the learning paiameters, and/oi use more hidden units The above pro 
ceduie is cilled met emenial backprop learning because the weights are updated 
after every pic sent atioii of an input pattern This incremental updating is more 
desnabk duo to two reasons (1) it requires less storage ,and (2) it makes the 
search path in Uic weight space stochastic (here, at each time step, the input 
vector X i& drawn at random), which allows for a wider exploration of the search 
space uid, potentully leads to better quality solution 


6 5 Backpropagation Enhancements and 
Variations 

In general, learning with backpropagation is slow Typically, this is due to the 
characteristics of the error surface[10] The surface is characterized by numerous 
flat uul Itcc-I. Kg.ons In addition, it has many tioughs that ate flat m the 
ditt ctioll ol scauli Some times when backpiopagation converges, it may converge 
to a local mimimim of the ciitonon function This fact is true in any gradient- 
descent based liaining rule when the surface rs being searched m nonconvex, 

1 c . li admits local minima The local minima problem cim be explained an 
shown in figure 6 8 Suppose one staited with a weight set for the network 
coirespondmg to point P Due to gradient descent learning, the mmimnm one 
can encounloi is the one at Mr not that at M, Mr rs called a tal m.mmum 
and corresponds to a parUal solution tor the network in response to the traim^ 
data M, IS the global minimum point which one have to ge , nn ess 
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Figure 6 8 Local Minimum 

die taken to escape from Mi, Mg will never be reached Many enhancements 
of and variations to backprop have been proposed These are mostly heuristic 
modifications with goals of incieased speed of convergence, avoidance of local 
minima, and/ oi impiovemeiit in the network’s ability to generalize Some of them 
aiG presented here which axe mainly considered by the author in his project 

6 5 1 Weight Initialization 

Owing to its gradient descent nature, backprop agation is very sensitive to initial 
conditions If the choice of the initial weight vector w° happens to be located 
ncaiei to the solution point, then the convergence of back prop is fast On the 
other hand, backprop converges very slowly if staxts the search in a relatively 
flat region of the ei ror surface In practice, the weights are normally initialized to 
small zero moan random values The motivation for starting from small weights 
IS that large weights tend to prematuiely saturate units in a network and render 
them insensitive to the learning process On the other hand, randomness is 


i 


I 
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mUodutcd as a, symmetay bredang mechamsm it prevents units from adopting 
siimlu [unctions and bccomming redundant [ 10 ] 

6 5 2 Learning Rate 

J ho convcigence speed of backpiop is dnectly related to the learning rate pa- 
1 xincioi rj If 17 IS small, the seaich path will closely approximate the gradient 
path, but conveigcnce will be veiy slow due to the large number of update steps 
needed to lench a minima On the other hand if t) is large convergence initially 
will be veiy fast, but the nlgoiithm will eventually oscillate and thus not reach a 
imiiiniLiiii In geneial, it is desnable to have large steps when the search point is 
1 11 away fioni a minimum, with decreasing step size as the search approaches a 
mini mum Many heuristics have been proposed so as to adapt the learning rate 
uitoin itu ally In [ 10 ] used a method where the learning rate for a given weight 
Wi IS to set to ar}t{t) if and have the same sign, with a > 1 if 

the ]) ntial dciivatives have different signs, then a learning rate of b 77 »(t) is used, 
with 0 < 6 < 1 A similai, theoietically justified method for increasing the con- 
vciguKc speed of inciemental gradient descent search is to set 7]{t) = 'q(t — 1) 
if VjE/(i) has the same sign as VE{t — 1), and r]{t) — r}{t — l)/2 otherwise 

6 5 3 Momentum 

Anothci simple approach to speed up backpropagation is through the addition of 
1 moiiK Ilium teim to the light hand side of the weight update rules m equations 
0 19 and 6 28 Ilcie, each weight change Aw* is given some momentum so that 
il accclcuitcs m the aveiagc downhill diiection instead of fluctuating with every 
change m the sign of the associated partial derivative The addition of 

momentum to gradient search is stated formally as 

- 1 ) ( 6 ^ 3 ) 

dw^it) 

where /3 is a momentum rate normally chosen between [0,1] and Atii,(t — 1) = 

'll ) — 1) Equation 6 33 is a special case of multistage gradient methods 
that have been proposed for acceleiatmg convergence and escaping local minima 
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I'lguu 6 9 Plots of f(nel) = tanh(net) and it’s derivative f’(n.et) 

6 5 4 The Flat Spot Problem 

Duuiig Ihc tiammg, if a unit in a multilayei network receives a weighted signal 
net with d larg( magnitue, this unit outputs a value close to one of the saturation 
levels of its activation function If the coiiesponding target value is substantially 
dilTeieni fioin that of the saturated unit, one can say that the unit is incorrectly 
satuiatcd oi has cntcied a flat spot When this happens, the size of the weight 
update due 1o backprop will be very small, even though the error is relatively 
laige, ind it will take an excessively long timefoi such incorrectly saturated units 
to itveisc Ihcu stxtcs Ihis situation can be explained by referring to figure 
6 9 wluie the activation function f{net) = tanh{l3net) and its derivative 
f'inet) = /3[1 - Pinet)] aic plotted foi /? = 1 Here, when ('net) has 
laigc iniguitudc, appioaches zeio Thus the weight change approaches zero 
in ccpiations 6 L9 and 6 28 even when there is difference between the actual and 
desired output lor a given unit A simple solution to the flat spot problem is 
to bias the derivative of the activation function, that is, replace and ff^ va. 
equations 6 20 and 6 27 by + e and /' + e respectively (a typical value for e 

IS 0 1) 



7 Conclusions 


Tlii'^ thesis wis focused on one of the problems concexniiig the trading of renew 
able cneigies, that is, on the fluctuations of the renewable production Those 
fluctuation pioblems apply specially to wind and photovoltaic power production 
A mnkctci who buys such ’’fluctuating lenewable energy”, could possible order 
ciicigy lioiii ( ojiV( iitional powei markets to compensate the fluctuation of the 
icnewabk piodiution To know how much he should order, he would need a 
good foiecast foi the renewable pioduction in the commmg hours/days Within 
the scope of tins work diffeient power foiecast models for wind and photovoltaic 
eneigy wcie established and compared 

Two wind powc i forecast models were made, a fuzzy based model and a multiple 
legiession base d model The results from the fuzzy model were acceptable With 
some advancement the model would probably give a better forecast More of the 
influencing f ictois could foi example be taking into account and a better fuzzy 
rule base could be built Ihe regression model gave considerable better results 
than I he fuzzy model But also here some advancement could be made In [21] 
it IS imiilioiicd, that the wind power has relationship to the cube of the wind 
speed flhis could be lakcn into account in the regression model in a hope for 
belter results And, like for the other model, some improvement could possibly 
be reached thiough higlier number of input variables 

For the photovoltaic power forecasting two models were constructed, a neural 
network based model and a multiple regression based model The results from 
both the models were quite good The neural network model gave slightly better 
results than the regression model The used data included recognizable noise and 
It would be interesting to test the models with more accurate data 
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111 ilu ‘1 IIk sis it lias be shown that with suitable models relatively good renewable 
powt 1 fo] cc ists can be made Such models could be extremely helpful for a power 
nuikctci, who IS concentiatcd on buying and selling renewable energy With the 
h( Ip ol such mode Is the lenewable eneigies could possibly be able to play moie 
impoilanl loh on the powei maikets m the future, than they have done until 


now 



A. Additional Figures and Tables 

A 1 Wind Power Forecasting 


S No 

Wind Speed 
m/s 

Actual Power 
kW 

Estimated Power 
kW 

Error 

1 

12 

48 6 

45 013560 

0 0"7379493 

2 

47 

63 5 

117 536520 

0 850968 

1 

27 

34 4 

36 226521 

0 0530965 

4 

55 

69 3 

156 288002 

1 255238 

5 

6 1 

104 9 

175 042783 

0 668663 

6 

46 

90 2 

113 456706 

0 257834 

7 

6 3 

58 5 

185 8826813 

2 177481 

8 

38 

52 7 

67 1450844 

0 2741002 

9 

24 

20 9 

35 595315 

0 7031251 

10 

49 

124 2 

142 250257 

0 14533 


lablc A 1 Wind power forecasting results using fuzzy model with regularity 
lufoimation ’’regular” 
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Piguie A 1 18th October Oster wind plant data 
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Appendiv A Additional Piguies and Tables 
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Figuie A 2 14tli November Oster wind plant data 
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Figuu A 3 14tli Novembei fuzzy based model result with regularity information 
’’little legular” 



Figiue A 4 1 4th November fuzz;y based model result with regularity information 

”veiy legular” 
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A 2 Photovoltaic Powgi* Forecasting 


V ui iblc 

6 

Sb 

Fb 

Correlation with Power 

R idi ition 

1 ( inp M iluic 
Wind Speed 

0 8788 

0 1184 

0 05608 

0 01077 

0 011781 

0 01329 

6651 304 
101 0244 
17 8029 

0 96943 

0 68491 

0 59891 


(T Ra) — ^ 

Rp(RaT) = 0 972415 

^p(WRaT) ~ 9 376E-004 

Fk2 = 59M937 


Table A 2 20tii June data multiple regression analysis 


Vunbh 

b 

Sb 

Fb 

Correlation with Power 

Radiation 

0 91547 

0 0444 

424 994 

0 81403 

luupai ituic 

0 12543 

0 03627 

11 9556 

0 6529 

Wind Speed 

0 00863 

0 0268 

0 10376 

0 26614 


j52 

■^^PiTRa) 


0 0066 Rp(RaT)= 0 8m 


^PjWRaT) 


1 963E-004 


Fr2 = 402 4229 


1 able A 3 9th January data multiple regression analysis 
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r iguie A 5 Solar Data for June 20th 1994 
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Time 



Figuie A 6 Solar Data for June 21st 1994 
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Figure A 7 Solar Data for June 22nd 1994 


Radialion(W/m ) 



Figuie A 8 Solar Data for January 9th 1994 
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Time 


Figuie A 9 Solar Data for January lOth 1994 
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Figure A 10 Solar Data for January lltb 1994 
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riguic All Phase I 9ih January data neural network results 



I iguie A 12 Phase I 20th june data neural network results 
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Figiuc A J 4 Phase I 22nd june 9th hour data neural network results 
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Figuie A 15 Phase II 22nd june 9th hour data neural network results 



I- iguie A 16 Phase I 9th January data multiple regression results 
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Figuic A 17 Pliase I 20ih june data multiple regression results 
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Figure A 18 Phase I 20th juiie 14th hour data multiple regression results 
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Figuie A 19 Phase I 22nd june 9th hour data multiple regression results 



F,gmo A 20 Pha.e 11 22nd june 9lh hour data multiple legression lesute 




Bibliography 


[1] Alexiadis, M C , Dokopoulos, P S (IEEE member), Sahsamanoglou, H 
S , Wind Speed and Potoer Forecasting based on Spatial Correlation Models, 
IFPE iiansaciioiis on Energy Conversion, Vol 14, No 3, September 1999, 
pp b36 842 

[2] Annan, R II , Her wing, L 0 Photovoltaic industry and electric utility roles 
in ike U S photovoltaic program, Proceedings of the Biennial Congress of 
the International Solar Energy Society, Denver, p 331, 1991 

[3] Applcbaum, J The quality of load matching in a direct coupling photovoltaic 
system, IEEE Trans on Energy Conversron, Vol 2, No 4, pp 534 541, 1987 

[4] Beyer , et a H G Short term prediction of wind speed and power output of a 
wind turbine with neural networks, Proed of the 5th European Wind Energy 
Conlcrence, EWEC94, Thessaloniki, Greece, 1994, pp 349 352 

[5] Bock, U , Nisscn, J Standardisiert Lastprofile fur Haushalte und 
Kkmqiweibc, Lngergiewirtschaftliche Tagesfragen, Jg 49 (1999), Heft 9, 

p606 610 

[6] Bossanyi, E A Stochastic wind prediction for wind turbine system con 
trol, Pioccl of the 7th BWEA Wind Energy Conference, 27 29 March 1985, 

pp 219 226 

[7] Chatfield, C Statistics foi Technology, A comce in Applied Statistics, 

[8] Edwards, A L Multiple Regression and the analysis of variance and co 
variance, W II rreeman & Co , pp 1-66, 1979 


100 



Page 101 


[9] I Is '1111 Wind powet piedtciizon tool in central dispatch centres^ Progress 
Rcpoit 3 ol Ihc LU project JOU2 CT92 0083, October 1993 

[10] 11 issoun, ]V1 II ] undarncHtal*^ oj A.itificial Meuicil Metwovhs^ Prentice Hall 
ol Indn piivite limited New Delhi, 1998 

[11] Ilaykiii, S Neural Nelwoiks, Pienlice Hall, Inc , New Jersey, 1999 

[12] Iliyama, T, , Kitabayashi, K Neural Network Based Estimation of Maxi 
mum Powei Geneiation from PV Module Using Environmental Information, 
IEEE Tians'ictions on Eneigy Conveision, Vol 12, No 3, September 1997 

[13] Jewell, W T. , Uniuh, T D limits on cloud induced fluctuation in pho 
tovoliau ginctaiwn, IEEE fiansactions on Energy Conversion, Vol 5, No 1, 
Maich 1990 

[14] Kaltschmitt, M , Wies(Hisg), A Erneuerbare Energien, System technik, 
Wiitsoliaftlichkeit imweltaspekte ISBN3 540 59362 4 Springer verlag Berlin 
Ileidelbeig, 1996 

[15] Kaimiotcvkis, G , Nogaiet, E , Staviakakis, G Wind power forecasting using 
advanced neural network models, IEEE Transactions on Energy Conversion 
Vol 11 No 4 pp 762 767, Decembei 1996 

[16] Kosko, B Newal Networks and Fuzzy Systems A Dynamic Systems Ap 
pioach io Machine Intelligence, Englewook Cliffs, NJ Prentice Hall, 1992 

[17] Olla Fifth renewable oidei foi England and Wales, Birmingham, Septem 
bei 1998 

[18] Ponce dc Lcao, M 3 , Matos, M A Fuzzy Models for Producers fiom 
Nalurla Resow ces, IEEE pioceding 1998 

[19] Ross, T J Fuzzy Logic with Engineering Applications , McGraw Hill 1995 

[20] Sakk, E , Thomas, R , Zimmeiman, R Power System Tournaments for a 
Deregulated Envnonmeni, System Sciences, 1997, Proceedings of the Thirti 
eth Hawaii International Conference on System Sciences, p681-686 



Page 102 


[2J] Saoud, 7 S , Jenkins, N , Simple wind farm dynamic model lEE Proc 

Gcnci alien Tiansimission Dislmbjitipp, VqI 142, No 5 September 1995 

I 

[22] Scliicirci, G J, Boots, M G, Martens, J W Voogt, M H Tradable 
Gntn Ceibificaies A new market based incentive scheme for renewable en 
cigy Iniioduction and analysis, Report from Netherlands Energy Research 
Foundation EGN, Maich 1999 

[23] Schwaitz, R J Photovoltaic Power Generation, Procedmgs of the IEEE, 
Vol 81, No 3, March 1993 

[24] Shrihui, Li, Don, C , Wunsch , O’Hair, E , Giesselraann, M G Neural 
Nctwoil foi Wind Power Geneiation with Compressing Function, 

[25] 1 widcll, J Renewable Enetgy Implementation and benefits TEE 2nd Inter 
national Conference on Advances in Power System Control, Operation and 
Management, December 1993, Hong Kong 

[26] Wilson, R E , Spera, D A Aerodynamic Behavior of Wind Turbines, in 
Wind Turbine Technology, AMSE Piess, chap 5, pp 215 282, 1995 

[27] Winter, C J , Sizmann, R L , Vant Hull, L L Solar Power Plants, Fun 
dainentals. Technology Systems and Economics Springer Verlag Berlin 
lleidelburg 1991 

[28] /nuada, J M Introduction to aitificwl newal systems StPaul[ua] West 
I'^ulilishmg, 1992 



^ y 

Ee/2^0efr»f 


D-\ 

^tl^ooolfn 


p^iioW 


T^35'P 


T^rJOc/M-ci-loV^ ^ 



130924 




130924 


Date^iBp 

This book Is to be returned on the 


date last stamped 



